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We observe a narrow enhancement near 2m v in the invariant mass spectrum of pp pairs 
from radiative J/tp — > -ypp decays. No similar structure is seen in J/ip —> Tv°pp decays. The 
results are based on an analysis of a 58 million event sample of J/ip decays accumulated 
with the BESII detector at the Beijing electron-positron collider. The enhancement can be 
fit with either an S- or P-wave Breit Wigner resonance function; in the case of the S-wave 
fit, the peak mass is below 2m p at M = 1859J^q (stat)jt,g (sys) MeV/c 2 and the total width 

is r < 30 MeV/c 2 at the 90 percent confidence level. These mass and width values are not 



consistent with the properties of any known particle. 

§1. Introduction 
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There is an accumulation of evidence for anomalous behavior in the proton- 

>^ ' antiproton (pp) system very near the M pp = 2m p mass threshold. The observed 

^ ■ cross section ^ for e + e~ — ► hadrons has a narrow dip-like structure at a center of 

Qh! mass energy of yfs ~ 2m p c 2 . The proton's time- like magnetic form- factor, deter- 



mined from high statistics measurements of the pp — ► e + e~ annihilation process, 
exhibits a very steep fall-off just above the pp mass threshold 2 >. The authors of 
ref. 1 attribute these features as being due to the effects of a narrow, S'-wave triplet 
pp resonance with J PC = 1 , a mass of 1870 MeV/c 2 , and a width of 10 MeV/c 2 . 
In studies of p annihilations at rest in deuterium, anomalies in the charged pion 
momentum spectrum from pd — > 7r~ir°p and 7r + 7r~n reactions 3 ) and the spectator 
proton spectrum from the pd — * 2ir + 3ir~p s process 4 ) have been interpreted as effects 
of a narrow sub-threshold resonance with properties similar to those of the proposed 
1 state listed above. There are no well established mesons that could be associ- 
ated with such a state. Belle 5 ) has reported observations of the decays B + — > K + pp 
and B° — ► D°pp. In both processes there are enhancements in the pp invariant mass 
distributions near Mpp ~ 2m p . 

The proximity in mass to 2m p is suggestive of nucleon-antinucleon (NN) bound 
states, an idea that has a long history. In 1949, Fermi and Yang 6 ) proposed that the 
pion was a tightly bound NN state. Although this turned out not to be correct, the 
train of thought started by this paper had enormous consequences. In 1956, after the 
discovery of the strange particles A and K° , Sakata 7 ) expanded this picture and used 
"fundamental baryon triplets" comprised of (p, n, A) and their antiparticles (p, n, A) 
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to make both pions and kaons. It was subsequently realized that the underlying 
mathematics of this model was that of the SU (3) unitary group 8 ) , a realization that 
inevitably led to the quark model. In 1961, Nambu and Jona-Lasinio 9 ) introduced 
a dynamical theory based on chiral invariance that also considers mesons as baryon- 
antibaryon composites. In this model, in addition to a low-mass pion, there is a 
scalar pp composite state with mass equal to 2m p . Here again, these ideas have been 
superseded by the quark model 10 ), and the scalar composite state is the famous a 
meson that is a main subject of discussion at this meeting. 

Although in both examples cited above the original motivation for nucleon- 
antinucleon composites is gone, the possibility of bound NN states with mass near 
2m p , i.e., p-p analogues of the deuteron and generally referred to as baryonium, 
continues to be considered u ). An investigation of low mass pp systems with different 
quantum numbers may help clarify the situation. 

In this talk I report on a study of the low mass pp pairs produced via radiative 
decays in a sample of 58 million J/ip events accumulated in the upgraded Beijing 
Spectrometer (BESII) located at the Beijing Electron-Positron Collider (BEPC) at 
the Beijing Institute of High Energy Physics. In this reaction, the pp pair is pro- 
duced in an environment that is free of any other hadrons. Moreover, charge-parity 
conservation insures that the pp system has C = +1. 



§2. Experimental considerations 



BESII is a large-solid-angle magnetic spectrometer that is described in detail 
elsewhere 12 ^. For this analysis we use events with a high energy gamma ray and 
two oppositely charged tracks that are identified as protons by their specific ioniza- 
tion (dE/dx) in the tracking chamber. Since antiprotons that stop in the material 
of the detection systems can produce annihilation products that are reconstructed 
elsewhere as 7 rays, no restrictions are placed on the total number of photons in the 
event. 

We subject the selected event candidates 

to four-constraint kinematic fits to the hy- 250 1 - 

potheses J/ip — > 'ypp and J/ip — ► r yK + K~. 
For events with more that one 7, we select 
the 7 that has the highest fit confidence level. 
We select events that have fit confidence level ° 
CL 1P p > 0.05 and reject events that have | ( 
CLjk+k- > CLypp. m 

Figure 1 shows the pp invariant mass dis- 
tribution for surviving events. The distribu- 
tion has a peak near M pp = 2.98 GeV/c 2 that 
is consistent in mass, width, and yield with ex- 
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Vc - ► PP 5 a broad 
2.2 GeV/c 2 , and a 
narrow, low-mass peak at the pp mass thresh- 
old that is the main subject of this talk 14 ). 



enhancement around M pp - ~ 



Fig. 1. The pp invariant mass distri- 
bution for the J/ip — ► 7pp-cnriched 
event sample 
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2.1. Backgrounds 

Backgrounds from processes involving charged particles that are not protons 
and antiprotons are negligibly small. In addition to being well separated from other 
charged particles by the dE/dx measurements and the kinematic fit, the protons and 
antiprotons from the low M pp - region tend to stop in the material in front of the elec- 
tromagnetic shower detector, where they have very characteristic responses: protons 
do not produce any signals in the shower detector while secondary particles from 
antiproton annihilation usually produce large signals. This asymmetric behavior is 
quite distinct from that for K + K~, 7r + 7r~ or e + e~ pairs, where the positive and 
negative tracks produce similar, non-zero responses. The observed shower counter 
energy distributions for the low-mass J/ip — > jpp events closely match expectations 
for protons and antiprotons and show no evidence for contamination from other 
particle species. 

There is, however, a large back- 
ground from J/ip —> 7r°pp events with an 
asymmetric 7r° — > 77 decay where one of 
the photons has most of the 7r° energy. 
This is studied using a sample of J/ip — > 
7T°pp decays reconstructed from the same 
data sample. For these, we select events 
with oppositely charged tracks that are 
identified as protons and with two or 
more photons, apply a four-constraint 
kinematic fit to the hypothesis J/ip — ► 
77P]5, and require CL 77P p > 0.005. For 
events with more than two 7's, we se- 
lect the 7 pair that produces the best fit. 
In the M 77 distribution of the selected 
events there is a very distinct ir° signal; 
we require |M 77 - M n o\ < 0.03 GeV/c 2 
(±2<r). The distribution of events vs. 
Mpp — 2m p near the M pp - = 2m p thresh- 
old, shown in Fig. 2(a), is reasonably well 

described by a function of the form /bkg(^) = N(5 1 / 2 + cqj 3 / 2 + a2<5 5 / 2 ), where 
6 = Mpp — 2m p and the shape parameters a\ and 02 are determined from a fit to 
simulated MC events that were generated uniformly in phase space. This is shown 
in the figure as a smooth curve. There is no indication of a narrow peak at low pp 
invariant masses. Monte Carlo simulations of other J/ip decay processes with final- 
state pp pairs indicate that backgrounds from processes other than J/ijj — > ir°pp are 
negligibly small. 

The Mpp — 2m p distribution for the 7r°pp phase-space MC events that pass the 
^ypp selection is shown in Fig. 2(b). There is no clustering at threshold; the smooth 
curve is the result of a fit to /bkg(^) with the same shape parameter values. 
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Fig. 2. The M P p — 2m p distribution for (a) 
selected J/tp — > ir°pp decays and (b) MC 
J/4> —> n°pp events that satisfy the ypp 
selection criteria. The smooth curves are 
the result of fits described in the text. 
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§3. Results 

Figure 3(a) shows the near-threshold M pp — 2m p distribution for the selected 
J ftp — > "ypp events. The solid curve shows the result of a fit using an acceptance- 
weighted S-wave Breit-Wigner (BW) function 15 ) to represent the low-mass enhance- 
ment plus /bkg(5) to represent the background. The mass and width of the BW signal 
function are allowed to vary and the shape parameters of fbkg(d) are fixed at the 
values derived from the fit to the ir°pp phase-space MC sample 16 **. This fit yields 
928 ± 57 events in the BW function with a peak mass of M = 1859 MeV/c 2 
and a full width of T = 0+^ MeV/c 2 . Here the errors are statistical only. The fit 
confidence level is 46.2% ( X 2 /d.o.f. = 56.3/56). 

Further evidence that the peak mass 
is below the 2m p threshold is provided in 
Fig. 3(b), which shows the M pp —2m p dis- 
tribution when the kinematic threshold 
behavior is removed by weighting each 
event by qo/q, where q is the proton mo- 
mentum in the pp restframe and qo is the 
value for M pp = 2 GeV/c 2 . The sharp 
and monotonic increase at threshold that 
is observed in this weighted histogram 
can only occur for an S'-wave BW func- 
tion when the peak mass is below 2m p . 

An S'-wave pp system with even C- 
parity would correspond to a CT + pseu- 
doscalar state. We also tried to fit 
the signal with a P-wave BW function, 
which would correspond to a ++ (qP) 
scalar state that occurs in some mod- 
els 9 ): 11 ). This fit yields a peak mass 
M = 1876.4 ± 0.9 MeVc 2 , which is very 
nearly equal to 2m p , and a very narrow 
total width: r = 4.6 ± 1.8 MeVc 2 (sta- 
tistical errors only). The fit quality, x 2 /d.o.f. = 59.0/56, is worse than that for 
the S-wave BW but still acceptable. A fit with a L>-wave BW fails badly with 
X 2 /d.o.f. = 1405/56. 

3.1. Can this be the effect of any known resonance? 

In addition we tried fits that use known particle resonances to represent the low- 
mass peak. There are two spin-zero resonances listed in the PDG tables in this mass 
region 17 ); the 7/(1760) with M„( 1760 ) = 1760±11 MeV/c 2 and r^ (1760) = 60±16 MeV, 
and the vr(1800) with M 7r(1800) = 1801 ± 13 MeV/c 2 and r„.( 1800 ) = 2 10 ± 15 MeV. A 
fit with /bk g and an acceptance-weighted S-wave BW function with mass and width 
fixed at the PDG values for the 7/(1760) produces x 2 /d.o.f. = 323.4/58. A fit using 
a BW with the 7r(1800) parameters is worse. 
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Fig. 3. (a) The near-threshold M P p — 2m p dis- 
tribution for the 7pp event sample. The 
solid curve is the result of the fit described 
in the text; the dashed curve shows the fit- 
ted background function, (b) The M P p — 
2m p distribution with events weighted by 
Qo/q- 
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3.2. Production angle distribution 

For both the scalar or pseudoscalar case, 
the polar angle of the photon, 7 , would be 
distributed according to 1 + cos 2 7 . Figure 4 
shows the background-subtracted, acceptance- 
corrected | cos 6* T | distribution for events with 
M pV < 1.9 GeV and |cos0 7 | < 0.8. Here we 
have subtracted the | cos w o | distribution from 
the 7T°pp data sample, normalized to the area 
of /bkg(<5) for M P p < 1.9 GeV/c 2 to account for 
background. The solid curve shows the result 
of a fit for 1 + cos 2 # 7 to the | cos# 7 | < 0.8 re- 
gion; the dashed line shows the result of a simi- 
lar fit to sin 2 6* 7 . Although the data are not pre- 
cise enough to establish a 1 + cos 2 # 7 behavior, 
the distribution is consistent with expectations 
for a radiative transition to a pseudoscalar or 
scalar meson 18 ) . 

3.3. Systematic Errors 

We evaluate systematic errors on the mass 
and width from changes observed in the fitted 
values for fits with different bin sizes, with background shape parameters left as free 
parameters, different shapes for the acceptance variation, and different resolutions. A 
study based on an ensemble of Monte Carlo experiments for sub-threshold resonances 
demonstrates that, in the presence of background, the farther the peak is below 
threshold, the less reliable is the determination of its mass. The MC studies also 
indicate that the mass determination of a below-threshold resonance can be biased. 
We include the range of differences between input and output values seen in the MC 
study in the systematic errors. 

For the mass, we determine a systematic error of MeVc 2 . For the total width, 
we determine a 90% CL upper limit of r < 30 MeV/c 2 , where the limit includes the 
systematic error. 

3.4. Branching Fraction 

Using a MC-determined acceptance of 23%, we determine a product of branching 
fractions B(J/Y> -► 7 X(1859))£(X(1859) -> pp) = (7.0±0.4(stat)lJ;|(syst)) x 10~ 5 , 
where the systematic error includes uncertainties in the acceptance (10%), the total 
number of J/ip decays in the data sample (5%), and the effects of changing the 
various inputs to the fit ( + ^2%)- 

§4. Summary 

In summary, we observe a strong, near-threshold enhancement in the pp invariant 
mass distribution in the radiative decay process J/ip — > jpp. No similar structure 
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Fig. 4. The background-subtracted, acceptance- 
corrected | cos 1 1 distribution for 
J/ip — » 7pp-enriched events with 
M P p < 1.9 GeV/c 2 . The solid 
curve is a fit to a 1 + cos 2 9~, shape 
for the region |cos# 7 | < 0.8; the 
dashed curve is the result of a fit 
to sin 2 8-y. 
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is seen in J/ip —> TT°pp decays. The structure has properties consistent with either 
a J PC = CT + or ++ quantum number assignment and cannot be attributed to the 
effects of any known meson resonance. If interpreted as a single resonance, its 
peak mass is below the Mpp = 2m p threshold at 1859_~\g(stat) jt^syst) MeV/c 2 and 
its width is r < 30 MeV/c 2 at the 90% CL. These parameters are quite similar to 
those of the 1 state proposed in ref. 1, which strongly suggests that these states 
may be related. 
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